Mosaic spin models with topological order 
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We study a class of two-dimensional spin models with the Kitaev-type couplings in mosaic struc- 
ture lattices to implement topological orders. We show that they are exactly solvable by reducing 
them to some free Majorana fermion models with gauge symmetries. The typical case with a 4-8-8 
close packing is investigated in detail to display the quantum phases with Abelian and non-Abelian 
anyons. Its topological properties characterized by Chern numbers are revealed through the edge 
modes of its spectrum. 
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Introduction- The phenomenon of emergence (such as a 
phase transition) in a condensed matter system is usually 
understood according to the Landau symmetry-breaking 
theory (LSBT) 1 . There also exists a new kind of or- 
der called "topological order"^* 2 -^^ which cannot be de- 
scribed in the frame of the LSBT (e.g., fractional quan- 
tum Hall effect). The study of topological order in theo- 
retical and experimental aspects has been an active area 
of research^^^^^i 2 .^^.^ . Since local pertur- 
bations hardly destroy the topological properties, such 
topologically ordered states show exciting potential to en- 
code and process quantum information robustly 2 . There- 
fore it is significant and challenging to find more exactly 
solvable models showing various topological orders. 

In this Rapid Communication, the Kitaev's honeycomb 
model 2 - is generalized to the general mosaic spin models 
with different two-dimensional Bravais lattices of com- 
plex unit cells. Then we study the 4-8-8 case in detail to 
reveal the general and special properties of mosaic spin 
models. 

Our mosaic spin models are constructed with the basic 
block shown in Fig. HJa), which is a vertex with three dif- 
ferent types of spin couplings along x- (black solid link) , 
y- (blue dotted link), and z- (red double link) directions, 
respectively. In spite of the lattice symmetry, numer- 
ous spin models can be built based on this basic block. 
However, taking translational symmetry and rotational 
symmetry as much as possible into account, we regard 
each basic block as the common vertex of three isogons 
with m, n2 and 713 edges, so there are only four kinds of 
mosaic spin models 16 illustrated in Fig. [3b)-flje), called 
m-n2-n 3 mosaic models. 

Obviously, the 6-6-6 mosaic model is just Kitaev's hon- 
eycomb model^. Here, we remark that for given n\, 112, 
and 713 , there exist some unequivalent kinds of plane ar- 
rangement of x links, y links, and z links, but we only 
illustrate one of them in Fig. Q] The general Hamiltonian 
of all mosaic spin models reads as 



H 



E 

u—x,y,z 



E , 

(j,k)es(u) 
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where S(u) is the set of links with it-direction couplings. 
Perturbation theory study and abelian anyons - First, 
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FIG. 1: (Color online) (a) Basic block for mosaic spin models, 
which consists of three branches with x- (black solid link), y- 
(blue dotted link), and z- (red double link) type couplings, 
(b) 6-6-6 mosaic model, i.e., Kitaev's honeycomb model, (c) 
3-12-12 mosaic model, (d) 4-6-12 mosaic model, (e) 4-8-8 
mosaic model, e vortices lie on squares while m vortices lie on 
octagons, (f) and (g) The possible nonconstant terms of the 
effective Hamiltonian are obtained by flipping four spin pairs 
around a octagon (f) and a quatrefoil (g). (h) Kitaev's toric 
code model is the effective model of the 4-8-8 mosaic model 
when \J Z \ > \ J X \, \J y \. 



we study the 4-8-8 mosaic model as a typical illustra- 
tion in detail. To sec its topological properties, we first 
analyze its low energy excitations when the system is ini- 
tially spontaneously polarized with the strong Ising in- 
teraction Ho = —J z J2 z unks a j a k- Tne ground energy of 
Hq is Eq = —NJ Z , where N is the number of z links. 
For larger J z in comparison with J x and J y , we regard 
the transverse part V = H — H as a perturbation and 
then prove that the obtained effective Hamiltonian H e f f 
just describes Kitaev's toric code model 2 , which supports 
many topological issues of the original mosaic spin model. 

The ground eigenstates of Hq are highly degenerate, 
where each two spins connected by a z link can be ei- 
ther |TT) or III)- The fusion projection^ T| can map 
the Ith aligned spin pair \m, m) l to an effective spin 
\m) l (to =t or J.), i.e., T\ \m,m) l = \m) v We use the 
fusion projection and the Green function formalism to 
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calculate the effective Hamiltonian H e ff — X^o-^e// 
= E + TtV[l + G (S ) + G (-Bo) ^G (Eb)]*HT + • • • 
where Go {E ) = (£o — Hq) . We first obtain the con- 
stant zeroth order one, the vanishing first order and third 
order ones. Here, each terms cr x a x or cr^a v k in V flips two 
spins, increasing the energy by 4J 2 . Up to the second or- 
der perturbation, one V flips two spins and the other V 
flips them back, giving H^f f = -N(J% + J%)/(4J Z ) as 
a constant. As shown in Figs. QJf) andQJg), we take 
two o~ x a x and two cr^af from four V around one octagon 
or one quatrefoil in a particular order. Taking all the 
2 x 4! = 48 possible cases into account, we obtain the 
fourth order effective Hamiltonian 



of some irreducible representations, and each irreducible 



H 



(4) _ 
eff ~ 
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where the constant term was dropped, O and Q repre- 
sent the octagon and quatrefoil in the two-dimensional 
(2D) lattice. Up to a unitary transformation for spin 
rotation a v — > er 2 , a z — > a x , a x — > <r y , the above Hamil- 
tonian represents the Kitaev's toric code mode£. Thus 
the above fusion projection constructs a new Bravais lat- 
tice illustrated in Fig. [2(h) with the effective spins laying 
on its links. Considering Kitacv model @ possesses rich 
topological features characterized by m and e anyons, 
we conclude that m particles live on octagons while e 
particles live on squares in our model with original spin 
representation. 

Majorana fermion mapping with TLi-gauge symmetry 
- The 4-8-8 mosaic model consists of four equivalent 
simple sublattices, and a unit cell [see the green rhombus 
tablet in Fig. [U(a)] contains each of four kinds of vortices 
referred to as 1, 2, 3, and 4. According to Kitae\*£, we 
use the Majorana fermion operators to represent Pauli 
operators as a x = ib x c, a y = ib v c, and o~ z = ib z c, where 
Majorana operators b x ,b v ,b z , and c satisfy a 2 = 1, a(3 = 
-(3a for a, (3 € {b x ,b y ,b z ,c} and a ^ (3. Then, the 
Hamiltonian ([T]) can be rewritten as H = Y]^ k \G jkCjCk, 
where the operator-valued coupling Gjk = iJ u Zjk (u — 
x, y, z) if (j, k) 6 S(u); Gjk = when (j, k) £ S(u). Here, 
a link (j, k) determines a type of coupling u = u(j,k). 
Due to the vanishing anticommutator of 6" and b% , we 
have Zjk = —Zkj for j ^ k. 

For each site, the above-mentioned Majorana operators 
act on a AD space, but the physical subspace is only 2D. 
Thus we need to invoke a gauge transformation of I2 
group to project the extended space into the physical 
subspace through the projection operator D = b x b v b z c: 
I i/j) belongs to the physical subspace if and only if D \ ip) = 
With this physical projection, some eigenstates of H 
can be found exactly because Gjk lays on the center of an 
Abelian algebra generated by Zjk with [Zjk, H ] = and 
[Zj k ,Z m i] = 0. Since (Z jk ) 2 = l,Z jk = ibjb^ generates 
a Z2 group and its eigenvalues are Zjk = ±1. Therefore, 
{Zjk, k) € S(u), u = x,y, z} generate the symmetry 
group Z2 ® Z 2 ® • • • <8> I2 of the model; the whole Hilbert 
space is then decomposed according to the direct sum 



sector is characterized by {Zjk\(j,k) G S(u), 



i.e., the directions shown in Figs. [2Ka)-[2fc). 

Obviously, in each irreducible representation space, 
we can reduce the Hamiltonian (fTJ) into a quadratic 
form, which represents an effective Hamiltonian of free 
fermions for a given vortex arrangement. To character- 
ize the vortex configuration, we introduce square and oc- 
tagon plaquette operators Wp = a\a\a^o-\ and Wp 8 ^ = 



< ] =- n z ^ 8) = 



n z ^ w 



where dp(A) and dp(8) represent the sets of boundary 
links of square and octagon plaquettes with label p; the 
(j, k) links are ordered clockwise around the plaquette. 
The operators WjP (j = 4, 8) commute with each other, 



0, W^ )2 = W^ 1 ' = I, and thus each pla- 
quette operator has two eigenvalues w p = ±1. A plaque- 



,(8)2 



tte with Wr, 



1 is a vortex-free plaquette while w p 



-1 



corresponds to a vortex. In the following we will show 
that different arrangements of vortices result in different 
phase graphs and different energy spectrums. 

4-8-8 mosaic model in different vortex-occupied sec- 
tors - Let us denote the site index j in detail by 
(s, A), where s refers to a unit cell, and A to a posi- 
tion type inside the cell. The Hamiltonian then reads 
H = J2 s ,\,t,a ( - ! sX,t t iC s xCt tJ ,/2. Due to the translational 
invariance of the lattice along the unit direction vec- 
tors ni = (1,0) , n 2 = (0,1), G S A,t M actually depends 
on A, /j and t — s, and thus exp [iq ■ (r t — r s )] G s \,t^ — 
exp (iq • r t ) GoA,t/j- To study the spectral structure of 
the system, we invoke the generic fermion operator a q M 



J\2N where N is the total number of the 
unit cells and Op^o^ + a^ x a p ^ = Spqd^x- Ga m (q) 
is the Fourier transformation of Go\,tfi- In the momen- 
tum space, the fermion representation of the Hamiltonian 
reads 



h = \Y. a \ g ^) a <- 



(4) 



Case I: In the vortex-free (VF) sector, we choose a 
particular direction (zjk = +1 or —1) for each link 
[see Fig. [UJa)], so that translational symmetry holds 
and Wp 4 ^ = Wp 8 ^ = 1 for all plaquettes. Since = 
(al 1, a q 2 , a q 3 , a q 4 ), we have the 4x4 spectral matrix 
G (q) = G VF or 

~ _ / J x a v -iJ y a x +iJ z a\ , , 

GvF - \ ij y o x - ij z0 t j x0 -y J ' (5) 

where a — diag.[exp(—iq2), — exp(igi)], q\ = q • nj, = 
q- n 2 . 

The single particle spectrum e (q) = —e (— q) is given 
by the eigenvalues of the spectral matrix G (q) . An im- 
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FIG. 2: (Color online) (a)-(c) 4-8-8 mosaic spin models in 
(a) vortex-free (VF) sector, (b) vortex-half occupied (VHO) 
sector, and (c) vortex-full occupied (VFO) sector, (d)-(f) The 
corresponding phase graphs of the above lattices with gapless 
phase B and gapped phases A: (d) VF, (e) VHO, and (f) 
VFO. 

portant property of the spectrum is whether it is gap- 
less, i.e., whether e (q) vanishes for some q. Obviously, 
the vanishing of determinant Det(GvF) enjoys the zero 
eigenvalues of Gvf- Then the gapless condition is 

J 2 + J 2 = J 2 . (6) 

As shown in Fig. Eld), the phase diagram of our model 
consists of three phases, the gapless phase B, which is ac- 
tually a conical surface, distinguishing from two gapped 
phases A z and A xy . Since the possible zero energy degen- 
erate points are (0, ±7r) and (±7T, 0) in the first Brillouin 
zone, we choose J x = J y = 1, J z = \[2, and qi = n to 
plot the profile graph of the energy spectrum with respect 
to q-2 G [— 7r, 7r] in Fig. G^b) by solid lines. The eigenval- 
ues of Gy f are chosen in the concourse {iv^cos (ga/4), 
±V2 sin (g>2/4)}. Thus in the vicinity of the energy de- 
generate points, the low-energy excited spectrum is ap- 
proximately linear. This property maybe helpful to study 
quantum state transfer problems^. 

Case II: We choose another particular direction for 
each link as shown in Fig. OJb), and the plaquettes with 

Wp 4 ^ = — 1 or Wp 8 ^ = — 1 are marked by blue shadings. 
In this vortex-half occupied (VHO) lattice, each unit cell 
contains eight kinds of sites, = (a q 1 , a q 2 , <z q 3 , a q 4, 

a q 5' a q 6' a q 7' a q s)> t ne corresponding 8x8 spectral ma- 
trix becomes 

/ J x cr y -U y Ux iJ z e^a f \ 

p iJyCTx Jx&y "iJzft 

UVH °- t iJ z (3 J x a y -iJ y a x 

\ -iJ z e lq2 a iJy(J x -J x a y / 

(7) 

where a = diag ^1, — e~ ,<?1 ^ , (3 = diag (e~ lqi , — 1^ , q[ = 

q ■ ni, q' 2 = q ■ n' 2 , n[ = (1, 1), and ri 2 = (-1, 1). The 
gapless condition for VHO lattice is 

Jx ^ Jy J !i Jy ^ J x Jz i Jz *~- Jx Jy (^) 



and the corresponding phase graph is plotted in Fig. [UJe). 
We notice that the same phase graph has been obtained 
by Pachosi for the Kitaev model. 

Case III: We choose the directions of links as shown in 
Fig. [2Jc) so that the translational symmetry still holds 
and w p 4 ^ = Wp = — 1 for every plaquette. The unit cell 
can be chosen as the same as the one in the VF sector, 
so do a, q%, and c?2- Therefore 

~ / J x a y -iJ y o x + ij z a\ 

GvFO ~ { iJ y a* - iJ^ ~J x ay ) ■ (9) 

The gapless condition is found as 

(Jx ~ Jzf < Jy < (Jx + Jz) 2 ■ (10) 

If J x , J y , J z > 0, we have J x < J y + J z , J y < J x + J z , 
Jz < Jx + Jy Thus in this case the phase diagram of 
our model is the same as that of Kitaev's honeycomb 
model. As shown in Fig. [2jf ) , the region within the red 
lines labeled by B is gapless. The other three gapped 
phases A x , A y , and A z are algebraically distinct. How- 
ever, the energy spectrum of the 4-8-8 mosaic model is 
more complex than that of the Kitaev model. When 
Jx = Jy = Jz = 1 and <Zi = — <Z2 = q, the eigenval- 
ues of the single fermion are chosen in the concourse 
{-l/2±cos(g/2 + 7r/4), l/2±cos(q/2-7r/4)}. Similarly, 
the different energy spectrums of mosaic spin models im- 
ply their different dynamic properties. 

In order to see the stability of the ground state in 
the VF sector, we compare the ground energy Eq — 
_ Eq £ q/ 2 of the VF lattice with that of the VHO 
and VFO lattice mentioned above. By choosing J x = 
J y = J z = 1, we find the ground energy per site is 
E y F = -0.80415, E ,vho = -0.75930, and E oyFO = 
-0.73631, so E G y F < E yHO < E yFO- The other 
cases can be studied similarly. Actually, as it was pointed 
out by Kitaev^, Lieb's theorem^ ensures that the VF lat- 
tice has the lowest energy to form a ground state. In the 
following, we will focus on the stable VF lattice and in- 
vestigate the nontrivial topological properties in the B 
phase. 

Topological properties of B phase in the pres- 
ence of magnetic field - The perturbation V = 
— (h x o-j + h y o V j + h z o-j) introduced by Kitaev^ can 
break the time-reversal symmetry. Then the nontrivial 
third-order term becomes H^A = k^2j k t (iZjiZ^i) CjCk, 
where k ~ h x h y h z / J 2 . As illustrated in Fig. EJa), 
the thin dashed arrows represent the effective second 
nearest-neighbor interactions between fermions induced 
by H e ff, and their directions denote the chosen gauge 
ZjiZki. When k = 0.025 , the changed profiled spec- 
trum is figured by dashed lines in Fig. [3jb). Therefore 
the system in the B phase acquires an energy gap in the 
presence of a magnetic field, which is helpful for protect- 
ing non-Abelian anyons. 

According to Kitaev^, the topological properties of a 
two-dimensional noninteracting fermion system with an 
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FIG. 3: (Color online) (a) Thin dashed arrows describe the ef- 
fective second nearest-neighbor interactions between fermions 
and the corresponding gauge induced by a magnetic field, 
(b) Profile graph of an energy spectrum with J x = J v = 1, 
J z = \pl along the q\ — it axis in the absence (solid lines) 
and presence (dashed lines) of a magnetic field, (c) Energy 
spectrum of the above system with finite size along the ri2 di- 
rection in the magnetic field. Two chiral edge modes crossing 
at E = correspond to Chern number ±1. 



energy gap are usually characterized by Chern number, 
which can be determined by observing the edge modes of 
the spectrum 2 -^. If the system illustrated in Fig. 
is finite along the ri2 direction while still periodic in the 
ri! direction, its energy spectrum is shown in Fig. [3Jc) 
with J x = J v = 1, J z = \/2, and k = 0.025. Then 
we can observe two edge modes (corresponding to two 
edges) crossing at E — 0. Since the Fermi energy lies in 
the central gap, only these two edge states around zero 
energy are relevant to Chern number—. We also notice 
that the two edge modes have a universal chiral feature 2 , 
i.e., even if the edges are changed, the energy curves of 
the edge modes do not change their tendencies, respec- 



tively. Therefore we conclude that the Chern number is 
±1 in the non-Abelian B phase. We also get zero Chern 
number in the Abelian phases A xy and A z with similar 
studies. Compared with Kitaev's honeycomb model and 
the 4-8-8 mosaic model with even cycles in the lattice, 
the 3-12-12 mosaic model with odd cycles spontaneously 
breaks time reversal symmetry to obtain Chern number 
±1 without applying a magnetic field 2 -^. 

Conclusion - We generalize Kitaev's honeycomb model 
to various mosaic spin models with translation and rota- 
tion symmetries and study the 4-8-8 case in detail. It is 
found that when \J Z \ 3> | J x \ , \ Jy \ , our model is equiva- 
lent to Kitaev's toric code model with Abelian anyons. 
Different vortex excitations result in different phase dia- 
grams with a gapless and gapped spectral structure. In 
the stable vortex-free case, the zero-energy Dirac points 
appear and the external magnetic field can induce an en- 
ergy gap. The nontrivial Chern number in B phase is 
obtained by studying the edge modes of the spectrum. 
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When this work is nearly finished, we notice that H. 
Yao and S. A. Kivelson have just studied the 3-12-12 
mosaic model in detail 2 ^. 
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